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Degenerate two-photon indirect absorption in silicon is an important limiting effect on the use of silicon 
structures for all-optical information processing at telecommunication wavelengths. We perform a full band- 
structure calculation to investigate two-photon indirect absorption in bulk silicon, using a pseudopotential 
description of the energy bands and an adiabatic bond charge model to describe phonon dispersion and 
polarization. Our results agree well with some recent experimental results. The transverse acoustic/optical 
phonon- assisted processes dominate. 

PACS numbers: 42.65.-k 



Silicon photonics^— has attracted widespread inter- 
est due to its potential application in ultracompact all- 
optical and optical-electronic devices, particularly for 
use at telecommunication wavelengths near 1.55 /xm. 
Many of these applications are based on the change in 
index of refraction with intensity, a third-order opti- 
cal nonlinearity. 4-6 Two-photon absorption, which can 
plague optical devices because of both immediate nonlin- 
ear loss and subsequent loss due to free-carrier absorp- 
tion, is also a third-order nonlinearity. However, in the 
interesting wavelength range of 1.2 /im to 1.7 /im the sum 
of energies from two photons is only sufficient to excite 
an electron across the indirect band gap, with the as- 
sistance of an absorbed or emitted phonon. So indirect 
degenerate two-photon absorption is of special interest 
for the development of optical devices r-~— and has been 
the subject of many investigations*^— 

Experimentally, the two-photon absorption coefficient 
j3 was first measured at 1.06 /im and 100 K by Rein- 
tjes and McGroddy^ and found to be 1.5 cm/GW. 
Thereafter a series of experiments at telecommunication 
wavelengths gave j3 in the range of 0.44 — 0.9 cm/GW. 6 
Bristow et aZ. 1 - measured /? as a function of pho- 
ton wavelength and found a strong dependence (see 
Fig. [I}. Anisotropy in the two-photon absorption was 
also investigated^ Theoretically, most studie d 18 ' 19 are 
based on the parabolic band approximation and phe- 
nomenological electron-phonon interactions. The exciton 
effect has also been investigated^ but only in quantum 
wells. Yet a full band structure calculation of the two- 
photon indirect absorption coefficient is still absent. 

In this paper we study the two-photon indirect ab- 
sorption spectrum in silicon. We perform a full band 
structure calculation using an empirical pseudopotential 
model (EPM) — for the electronic states, and an adia- 
batic bond charge model (ABCM)^ 2 - for the phonons, and 
present the anisotropy and temperature dependence of 
the absorption. The calculated two-photon absorption 
coefficient is in good agreement with some recent exper- 



imental results. 

For an incident laser beam with electric field E(t) = 
E^e" lut + c.c, the rate of the carrier density injection is 
written as 



(1) 



The superscripts indicate Cartesian coordinates, and re- 
peated superscripts are to be summed over; £ abcd = 
J2\± £x± d are the injection coefficients, with stand- 
ing for the contribution from the A'^-branch phonon 
emission (+) or absorption (— ) process. From Fermi's 
Golden Rule we have 
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Here k c and k v are the electron and hole wavevectors 
respectively; {a <H> b} stands for the same term as in 
the previous {}, but with the interchange a «-> b. The 
subscript c (v) is a band index for a conduction (va- 
lence) state; n, m identify intermediate states; is 
the electron energy at band n and wavevector fc, and 
fi^nrnk = Snk-^mk- The phonon energy is given by hQ qX 
for wavevector q and branch A (transverse acoustic (TA) 
and optical (TO), and longitudinal acoustic (LA) and op- 
tical (LO) branches), and N q \ is the equilibrium phonon 
number. The velocity matrix elements are v nm \~ = 



2 
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-g^-|mfc), with ff e standing for 

tron Hamiltonian, and the M n k a mk v x 
k v )\mk v ) are matrix elements of the electron-phonon in- 
teraction; H ep = J2x q ( ( l)( a <i*+ a -q\) i s tne electron- 
phonon interaction, with a q \ standing for the phonon 
annihilation operator. 

Due to the Oh symmetry in bulk silicon, the fourth- 
order tensor £ afccrf (with £ abcd = £, bacd and £ abcd = £ abdc ) 
only has three independent nonzero components: 
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the £f^ d share the same symmetry properties. 

For a quantitative calculation, we employ the EPM 
for the electronic states and the ABCM for the phonon 
states; all parameters are those used in the calculation of 
one-photon indirect absorption, 2 - where good agreement 
with experiment has been achieved. Within the pseu- 
dopotential scheme we determine the electron-phonon in- 
teraction, and then evaluate the matrix elements H^ p (q) 
using the calculated electron and phonon wavefunctions. 
At zero temperature, the calculated indirect band gap 
is Ei g = 1.17 eV, but it decreases with increasing 
temperature^ due to the electron-phonon interaction^ 
Thus we present our results as a function of the excess 
photon energy 2Tiu — E ig . In the numerical calculation, 
we use an improved adaptive linear analytic tetrahedral 
integration method to perform the six-fold integration, 23 
and choose the lowest 30 bands as intermediate states 
to ensure convergence. According to investigations on 
one— and two-photon 2 * indirect gap absorption within 
the parabolic band approximation, the neglect of the ex- 
citonic effects does not change the lineshape at energies 
more than a few binding energies above the band gap, 
and we neglect them here. However, since first principle 
studies of one-photon direct-gap absorption^ show that 
the electron-hole interaction plays an important role even 
for higher photon energies, further work is in order to in- 
vestigate its role in indirect gap absorption. 

We first compare our calculated indirect two-photon 
absorption coefficient with experiments listed in Fig. Q] 
for incident light linearly polarized along the [001] direc- 



tion; /3 



2/^ 



Here tir is the refractive index, 



j2n R ce ) 

c is the speed of light, and eo is the vacuum permittivity. 
We restrict our calculations to photon energies below the 
onset of one-photon indirect absorption; beyond this en- 
ergy divergences in Eq. ^ become problematic and a 
more sophisticated theory taking into account real popu- 
lation in intermediate states is necessary. Our results are 
consistent with the measured values of Bristow et alr^ 
The calculated j3 at 1.55 //m is around 1.0 cm/GW. 

When the laser is incident along a direction other than 
[001], or is not linearly polarized, other £ abcd contribute 
to the absorption and lead to its anisotropy Fig. [5] (a) 
shows the spectra of £ xxxx : ^ xx vv^ anc { ^ x v x v a t tempera- 
ture 4 K (red thin curves) and 300 K (black thick curves). 
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FIG. 1. (color online). Spectrum of at 300 K. The circles,^ 
uptriangles,— down triangles,— diamonds,— and squares^ 
are experimental results. The solid curve is our calculation. 
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FIG. 2. (color online). Detailed calculation results, (a) 
Spectra of £ xxxx , £ xxyy , and £ xyxy at 4 K (red thin curves) 
and 300 K (black thick curves), (b) Temperature depen- 
dence of {(T)/£(T = 0) at different photon energies Tiu)i: 
2hu)i—Eig — 0.1, 0.43, and 1 eV for ui, LU2 and LJ3 respectively; 
at K ui2 corresponds to a wavelength of 1.55 /^m. (c)-(d) 
Phonon-resolved spectra of £ xxyv at 4 K (c) and 300 K (d); 
black thick (red thin) curves are for phonon emission (absorp- 
tion) processes. 



When the excess photon energy is less than 1.2 eV, all 
injection coefficients increase with both photon energy 
and temperature, except that £ xx vy increases to a max- 
imum value and then decreases slightly with photon en- 
ergy at 300 K. Compared to the one-photon indirect in- 
jection coefficient, which is nearly quadratic in photon 
energy^ the two-photon indirect injection coefficients 
increase more slowly due to the explicit dependence on 
and the complicated frequency dependence of the 
transition matrix elements appearing in Eq. ([3]). 

To better understand the spectra, in Fig. [3] we also 
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plot the phonon-resolved spectra of at 4 K (c) 

and 300 K (d) as black thick curves. The LA phonon 
assisted process is least important at all temperatures, 
and gives a contribution less than 15%; the TA and TO 
phonon assisted processes give the most important con- 
tributions, and dominate at low and high photon ener- 
gies respectively. The crossover photon energy, where 
the TO phonon assisted process first dominates, moves 
to higher photon energy with increasing temperature be- 
cause the TA phonon has the smaller energy, and the 
TA phonon assisted process is more sensitive to temper- 
ature. At 300 K, the spectra of £ xxxx are given as red thin 
curves in Fig. [2](d). Only the TA phonon absorption as- 
sisted process gives a significant contribution. Therefore, 
the photon energy dependence of each absorption pro- 
cess is strongly determined by the properties of the me- 
diated phonons, and selection rules taking into account 
the phonon properties are necessary. 

In Fig.[5](b) we give the temperature dependence of the 
injection coefficient at different photon energies. When 
the temperature is higher than 200 K, the injection co- 
efficients depend approximately linearly on temperature, 
with the slope decreasing with photon energy. As with 
one-photon indirect injection^ the only temperature de- 
pendence comes from the equilibrium phonon number 
N q \, which is approximately linear with temperature at 
high temperatures. 

In contrast to the two-photon direct absorption coef- 
ficient, which has a cusp 30 i 31 that identifies the onset of 
absorption from the spin-split off band, the spectrum of 
fj in indirect absorption is very smooth. This is because 
the energy dependence at the onset of indirect absorp- 
tion is oc (hu) — Ei g ) 2 , instead of oc (hui — Eg) 1 / 2 for di- 
rect absorption. In contrast to one-photon indirect gap 
absorption^ the TA phonon assisted process here is at 
least as important as TO phonon assisted process. 

In conclusion, we have performed a full band calcu- 
lation of two-photon indirect absorption in bulk silicon. 
We find the two-photon indirect absorption coefficient is 
about 1 cm/GW at 1.55 and 300 K. The TA and TO 
assisted processes give the most important contributions 
at all temperatures. The injection coefficients depend 
linearly on temperature over 200 K. 
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